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Abstract. Stars form from dense molecular cores, and the mass function of these cores (the CMF) is often found to be similar 
to the form of the stellar initial mass function (IMF). This suggests that the form of the IMF is the result of the form of the CMF. 
However, most stars are thought to form in binary and multiple systems, therefore the relationship between the IMF and the 
CMF cannot be trivial. We test two star formation scenarios - one in which all stars form as binary or triple systems, and one in 
which low-mass stars form in a predominantly single mode. We show that from a log-normal CMF, similar to those observed, 
and expected on theoretical grounds, the model in which all stars form as multiples gives a better fit to the IMF. 
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1. Introduction 

^ The origin of the stellar initial mass function (IMF) is one of 
. the outstanding unsolved problems in astrophysics. As stars 
' form in dense molecular cores (see e.g. Ward-Thompson et al. 
. 1994; Kirk et al. 2005; Ward-Thompson et al. 2007), it might 
' well be expected that the IMF is related to the mass function 
of those cores (the CMF). This idea is supported by observa- 
tions of prestellar cores, which show that their mass functions 
are often similar to the IMF of Galactic field stars (Motte et al. 
1998; Testi & Sargent 1998; Johnstone et al. 2000; Johnstone 
et al. 2001; Motte et al. 2001; Johnstone & Bally 2006; Alves 
et al. 2007; Young et al. 2006; Nutter & Ward-Thompson 2007; 
Simpson et al. 2007). Further support is given by the observa- 
tion that Taurus may have both an unusual CMF (Onishi et al. 
2002) and an unusual IMF (Luhman 2004; see also Goodwin 
et al. 2004c), although Kroupa et al. (2003) show that the IMF 
in Taurus may be compatible with the field IMF. 

However, the relationship between the CMF and the IMF 
cannot be simple, as many, if not the vast majority, of stars form 
in binaries or higher-order multiple systems (see Goodwin & 
Kroupa 2005; Duchene et al. 2007 and Goodwin et al. 2007 
and references therein; see also Clark et al. 2007). Observations 
suggest that the binary frequency amongst young stars is 
higher than in the field (see Goodwin et al. 2007 and refer- 
ences therein) implying that binaries are destroyed by dynam- 
ical interactions in clusters (see Kroupa 1995a,b). However, 
Lada (2006) has argued that most M-dwarfs form as single 
stars, since the field M-dwarf binary fraction is relatively low 
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and there is no need to invoke dynamical destruction of low- 
mass binaries to form these (single) stars. The opposing view 
is argued by Goodwin & Kroupa (2005) and Goodwin & 
Whitworth (2007). 

If stars (or at least relatively high-mass stars) usually form 
in small-A^ multiples then there cannot be a trivial one-to-one 
relationship between the IMF to the CMF. Firstly, the mass of a 
core is distributed between a number of stars. Secondly, some 
stars are expected to be ejected at an early age from small-A^ 
multiples (e.g. Reipurth & Clarke 2001; Goodwin et al. 2007 
and references therein; see also Section 3). Thirdly, many bi- 
nary systems are expected to be destroyed in clusters (Kroupa 
1995a,b; Kroupa et al. 2003; Goodwin & Whitworth 2007; also 
see Goodwin et al. 2007 and references therein). Thus the CMF 
should relate most closely to the initial system mass function 
which, in turn, is modified by dynamical effects to produce a 
mixture of single and multiple systems. 

In this paper we examine the relationship between the IMF 
and the CMF, in particular we use the new results for the CMF 
in Orion from Nutter & Ward-Thompson (2007). In Section 2 
we review observations of the CMF, in Section 3 we present 
our general method, and in Section 4 we compare the IMFs we 
produce with the observations. 

2. Observations of the CMF 

The first observational link between the IMF and the CMF was 
made by Motte et al. (1998) in a millimetre study of the p- 
Ophiuchi molecular cloud. They found that the high-mass slope 
of the CMF matched that of the IMF. This result has been con- 
firmed for Ophiuchus (Johnstone et al. 2000; Young et al. 2006; 
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Simpson et al. 2007) and a number of other nearby clouds, 
including Orion (Motte et al. 2001; Johnstone et al. 2001; 
Johnstone & Bally 2006; Nutter & Ward-Thompson 2007), the 
Pipe Nebula (Alves et al. 2007), and Taurus (Onishi et al. 2002; 
however see Goodwin et al. 2004c), as well as for more dis- 
tant massive star-forming regions such as NGC 7538 and M17 
(Reid & Wilson 2006a,b). 

While the slope of the CMF seems to be consistent from 
region to region, the position of the peak of the CMF ap- 
pears to shift from ~ 0.1 Mq in nearby low-mass regions such 
as p-Ophiuchus (e.g. Motte et al. 1998), to a higher mass of 
~ 1 Mq in more distant and massive star-forming regions such 
as Orion (e.g. Nutter & Ward-Thompson 2007). Very mas- 
sive star-forming regions such as M17 show a flattening of 
the CMF at an even higher mass of ~ 8 Mq (Reid & Wilson 
2006a,b), though the data are incomplete before a turn-over 
is seen. Whether this is an intrinsic effect where the mass of 
the peak in the CMF is related to the mass of the stars being 
formed, or an observational effect caused by the blending of 
multiple sources at larger distances, is not yet known. 

3. The relationship between core and stellar mass 
functions 

We assume that the star formation properties of a core may be 
described by three basic parameters: the mass of the core Mc, 
the efficiency with which the core turns gas into star(s) e (so 
that the total mass of stars is eMc), and the number of stars 
formed within a core (the choices of A^, are discussed in 
section 4). We note that e and A^, may well be functions of Mc- 

There are then two basic distribution functions of the stars. 
The multiple system mass function (MSMF) is the mass distri- 
bution of the multiple systems produced by cores. The single 
star mass function (SSMF) is the mass distribution of all of the 
individual stars formed in all of the single and multiple sys- 
tems in all of the cores. It is the SSMF that will correspond to 
the (binary corrected) initial mass function (IMF; see below). 

Note that the MSMF will evolve due to several effects. 
Firstly, unstable high-order multiple systems may decay, pref- 
erentially ejecting single, low-mass stars (e.g. Reipurth & 
Clarke 2001; Bate et al. 2002, 2003; Sterzik & Durisen 
2003; Goodwin et al. 2004a,b; Delgado Donate et al. 
2004a,b; Hubber & Whitworth 2005; Umbriet et al. 2005). 
Secondly, binaries may be 'destroyed' by rather violent close 
binary-binary /binary-single interactions (Kroupa 1995a,b), and 
thirdly, wide, low-mass binaries may be 'disrupted' by more 
gentle impulses from passing stars (Goodwin & Whitworth 
2007). 

3.1. From a CMF to an IMF 

The procedure for generating an IMF (via a MSMF and SSMF) 
from a CMF is very simple. We randomly sample a core mass 
from a CMF (see Section 3.3 for our choice of CMF). This core 
then produces N* stars of total mass eMc (this is similar to the 
approach of Sterzik et al. 2001, however we do not constrain 
the IMF of the stars in any way except through the CMF of the 
cores). 



The masses of the A^, components in a multiple system are 
chosen randomly. In a binary system (A^, = 2), the masses of 
the primary and secondary are selected from a ffat mass ratio 
distribution (ie. one random number U[0, 1]). In a higher-order 
multiple the masses are distributed randomly, ie. A^, random 
numbers U[0, 1] are chosen and then the sum is normalised to 
unity to provide the mass distribution. 

The SFE, e, is chosen to provide the best fit to the canonical 
IMF (see below) and is assumed to be constant for all cores. It 
might be thought that SFE should depend on the mass of stars 
formed, as feedback energy increases with increasing stellar 
mass. However, the potential well from which gas must be re- 
moved by feedback also increases with increasing stellar and 
gas mass and so possibly the SFE is constant, or even increas- 
ing with mass. Given the uncertainties involved we make the 
simplest assumption possible that the SFE is constant. As will 
be seen, a good fit to the IMF can be obtained while making this 
assumption and there appears no need to appeal to a variation 
of SFE between cores of different masses. 

3.2. The canonical observed IMF 

We assume that the actual underlying IMF of stars has the 
canonical form (Kroupa 2002; see also Kroupa 2007) 
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with ai - 0.3, 0-2 = 1-3 and - 2.3 as the slopes, and Mo - 
O.O2M0, Ml = O.O8M0, Ml = O.5M0 and M3 = IOMq as die 
masses of the limits and turning points of the IMF. This form of 
the IMF matches well other recent determinations of the IMF 
(e.g. Chabrier 2003). 

The canonical IMF is corrected for the presence of unre- 
solved binary systems and therefore the IMF should follow the 
SSMF 

3.3. Forms of the CMF 

Our standard CMF is that determined for Orion by Nutter & 
Ward-Thompson (2007). We model this as a log-normal with 
dispersion criogj^M - 0.55 and a mean of yUiog,„M - 0.05 as 
illustrated in Fig.[T] This CMF is not too dissimilar to an IMF- 
like distribution with the slopes a remaining the same, but turn- 
over masses of ~ 8 times those in the canonical IMF (eqn.[TJ. 
We also note that this is very similar to the CMF of the Pipe 
dark cloud (Alves et al. 2007). 

Our fit to the Nutter & Ward-Thompson (2007) CMF is 
above the final two points which are below the completeness 
limit. We assume that the entire core mass distribution is mod- 
elled by a log-normal which would be expected on theoretical 
grounds (Padoan & Nordlund 2002,2004; Klessen & Burkert 
2000; Klessen 2001; Li et al. 2003; Jappsen et al. 2005) even 
in a wide variety of physical conditions (see esp. Padoan & 
Nordlund 2002; Jappsen et al. 2005). We also note that Padoan 
& Nordlund (2002,2004) tend to find a very steep decline in the 
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Fig. 1. The CMF of Nutter & Ward-Thompson (2007) with Va^ error-bars, the completeness limit of the observations is marked 
by the vertical dashed line. The data are fitted by a log-normal distribution with mean f^\ogi„M - 0.05 and criogj^M = 0.55 (dashed- 
dot line). Also plotted for comparison is the canonical IMF (eqn.[T]) with the same slopes {ai = 0.3, ai = 1.3 and ai - 2.3), but 
turn-over masses increased by a factor of 8 (Mo = O.I6M0, Mi = O.64M0, M2 - 4Mq and M3 = 8OM0). 



CMF below the peak which would be compatible with these 
observations. 

Clark et al. (2007) recently noted that the free-fall times 
of clumps of different masses are different and that low-mass 
cores collapse significantly faster than higher-mass cores. Thus 
an observed CMF must be constantly replenished with low- 
mass cores in order to retain a constant form. We assume that 
our CMF represents a 'snapshot' of the CMF. If star forma- 
tion occurs rapidly in clusters (in < 1 Myr, e.g. Elmegreen 
2000) then the observed CMF should represent the total CMF 
for masses above the peak of the CMF as the free-fall time for 
cores > IMq is a significant fraction of the cluster formation 
timescale (Clark et al. 2007, see their fig. 1). 

4. Results 

4.1. The fully multiple model 

Goodwin & Kroupa (2005) suggested that the observed prop- 
erties of multiple systems could be reproduced if each core 
produces 2 or 3 stars. Single field stars are then produced by 
the dynamical decay and destruction of multiple systems in 
young clusters (Kroupa 1995a,b; Goodwin & Kroupa 2005; 
Goodwin & Whitworth 2007; Goodwin et al. 2007 and refer- 
ences therein). 

To model the fully multiple scenario we assume that cores 
of mass eMc < Q.75Mq form entirely binary systems, and 
cores with eMc > Q.75Mq form multiple systems with a ra- 
tio of 3 : 1 binaries-to-triples. The SFE is chosen to give a good 
fit to the canonical IMF with e = 0.27. 



In this scenario the initial multiplicity fraction is unity. 
Single stars and brown dwarfs are produced by the destruction 
of many (especially low-mass) initially multiple systems (see 
Section 3). 

The result of the fully multiple model are illustrated in 
Fig. |2] This model produces a good fit to the canonical IMF 
for all masses except the very highest. The mass functions dip 
below the canonical slope of - 1 .3 at high masses due to the 
steep decline of the log-normal CMF at high massed 

The fully multiple model requires the dynamical destruc- 
tion (see e.g. Kroupa 1995b; Kroupa et al. 2003; Goodwin & 
Kroupa 2005; Goodwin et al. 2007) of significant numbers of 
low-mass binary systems in young clusters in order to change 
the initial binary fraction of unity to the field value. 

We note that in this model brown dwarfs are not primarily 
produced as single objects in cores ('star-like' formation, e.g. 
Padoan & Nordlund 2004), nor as ejected embryos from high- 
mass cores (the ejection hypothesis, e.g. Reipurth & Clarke 
2001). Instead they mainly form as the distant companions to 
M-dwarfs which are then disrupted. This is the scenario pro- 
posed by Goodwin & Whitworth (2007) as a major mode of 
brown dwarf formation. We note that this might be consistent 
with the idea that brown dwarfs form as a separate population 
of objects, possibly with a discontinuous IMF (Kroupa et al. 
2003; Thies & Kroupa 2007; Kumar & Schmeja 2007). 



' It might be expected that the actual CMF continues as a power- 
law decline rather than being fitted by a log-normal at high masses 
(see e.g. Padoan & Nordlund 2002,2004). 
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Fig. 2, The fully multiple model with Ni, - 2 for systems with eMc < 0.75 Mq and A^* =3:1 binary-to-triple ratio for 
eMc > O.75M0, and e = 0.27. The initial MSMF (dashed line), and the SSMF (red dotted line) compai-ed to the canonical IMF 
(thin solid lines). 



4.2. The Low-mass single star model 

Lada (2006) suggests that M-dwarfs tend to form as single 
stars, since most M-dwarfs in the field (roughly 55% by total 
number) are single stars. In this picture, destructive dynamical 
processes are unimportant for low-mass stars. However, Solar- 
type stars (and higher-masses) must still usually form in mul- 
tiple systems to fit the observed high multiplicity fraction of T 
Tauri stars (e.g. Mathieu 1994; Patience et al. 2002; Duchene 
et al. 2007 and Goodwin et al. 2007 and references therein). 

We model the situation where most low-mass stars form 
as singles (the low-mass single star model) by making Nt a 
strong function of the core mass. Cores with eMc < O.5M0 
form stars in a 2 : 1 single-to-binary ratio. Cores with masses 
0.5 < eMc/Mq < IMq have a 1 : 1 ratio of singles-to-binaries. 
And cores with eMc > IMq have a 3 : 1 binary-to-triple ratio 
(as in the fully multiple model). These probabilities roughly re- 
flect a combination of the low-mass field (M-dwarf, see Fischer 
& Marcey 1992, also Lada 2006) and intermediate-mass PMS 
binary fractions (see Duchene et al. 2007 and Goodwin et al. 
2007 and references therein). 

The results from this model are shown in Fig.[3]with a SFE 
of e = 0.15 which provides the best (but still not a good) fit 
to the canonical IMF. This model cannot be made to fit the 
canonical IMF well with any choice of e. 

The main problem with the low-mass single star model is 
that it provides no mechanism for forming brown dwarfs other 
than as single brown dwarfs in very low mass cores. As the 
CMF peaks at ~ IM0 there are very few cores below O.IMq 
and so the SFE must be made very low in order to shift the peak 
of the CMF to the peak of the IMF. This is avoided in the fully 
multiple model by producing brown dwarfs as companions to 



low-mass stars. Increasing the SFE to greater than 15% results 
in far too few VLMOs and also shifts the peak of the IMF to 
too-high a mass. We also note that a SFE of 15% suggests that 
a significant number of cores produce planetary-mass objects 
(which are not shown in Fig.O. 

There are two ways we might possibly escape from the 
problem of too-few brown dwarfs. Firstly, we may postulate 
a very large population of low-mass cores from which brown 
dwarfs can form (e.g. Padoan & Nordlund 2004). In this sit- 
uation the vast majority of cores would be below the peak of 
the CMF. However, the simulations of Padoan & Nordlund 
(2002,2004) suggest that the CMF should drop sharply below 
the peak, far faster than a log-normal distribution. Thus, there 
appears to be no good reason for expecting a large popula- 
tion of low-mass cores (see Goodwin & Whitworth 2007 for a 
number of other reasons why significant numbers of low-mass 
brown dwarf-forming cores would not be expected). 

Secondly, we may form brown dwarfs by ejecting low-mass 
embryos from massive cores, thus creating brown dwarfs (e.g. 
Reipurth & Clark 2001). However, large numbers of ejections 
per high-mass core would be required to create almost all of 
the brown dwarf population in this way, and large numbers of 
such ejections have consequences that are difficult to reconcile 
with observations (see e.g. Goodwin & Kroupa 2005; Goodwin 
et al. 2007 and Whitworth et al. 2007 and references therein; 
Goodwin & Whitworth 2007). 

The slope of the high-mass end of the IMF is also far too 
steep (~ 1.6 rather than the observed ~ 1.3). This is due to 
the change in the modes of fragmentation at eMc - 0.5Mq 
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and 1M(tQ. This is because at high masses the mass of stars is 
divided between 2 or 3 stars, while at lower masses it is divided 
only between 1 or 2. Thus altering the modes of fragmentation 
changes the slope from the IMF-like slope of the CMF to a 
steeper slope. The upper-mass slope only matches the slope 
of the CMF if fragmentation is independent of mass for cores 
above the knee in the IMF at O.SMq. 

The problem with the upper-mass slope can be alleviated 
somewhat by assuming that all stars form with the field binary 
fraction (rising from 33% for M-dwarfs to 60% for G-dwarfs). 
However, this solution conflicts with observations that the ini- 
tial binary fraction for stars > IMq is consistent with unity 
(Goodwin & Kroupa 2005; Duchene et al. 2007; Goodwin et 
al. 2007). These observations suggest that there must be a fairly 
rapid transition between 0.5 and IMq from a low to a high pri- 
mordial binary fraction which will result in too-steep an upper- 
mass slope of the IMF. 

We have assumed that the mass ratio of binaries is a flat dis- 
tribution. Biasing the mass ratio distribution to low-^ (ie. highly 
unequal mass systems) improves the problems at the high-mass 
end of the IMF slightly. If most high-mass cores produce one 
large star and one or two very low-mass stars, then the IMF at 
the high-mass end becomes more similar to the MSMF (as this 
is dominated by one of the stars). However, the mass ratio dis- 
tribution needs to be very biased for this to have a significant 
effect. 

The too-steep slope of the upper-end of the IMF can also 
be solved by assuming that the SFE increases with increasing 
core mass (in just the right way). However, we feel this solu- 
tion is unlikely as the SFE would have to be fine-tuned to give 
the correct slope and it would seem peculiar to postulate that 
low-mass cores produce stars at very low efficiencies (~ 10%), 
whilst higher-mass cores are able to convert more of their gas 
(~ 30%) into stars (the opposite of what might be expected 
from arguments based on feedback). 

5. Discussion 

We have shown that the observed mass functions of cores in 
Orion B (Nutter & Ward-Thompson 2007) can give rise to the 
IMF of stars. In particular we have shown that, to produce the 
stellar and sub-stellar IMF, the majority of these cores must 
fragment into multiple systems. However, there are a number 
of issues about cores and the CMF that are worth discussing in 
this context. 

It should be noted that it may not be fragmentation into 
'cores' in clusters that sets the IMF of stars. If competitive ac- 
cretion (see Bonnell et al. 2007 and references therein) is the 
dominant process, then the CMF at best acts to set the initial 
masses upon which competitive accretion begins to work. In 
such a scenario there would be little or no relationship between 
the CMF and the IMF. 

However, we would argue that the form of the CMF in dif- 
fuse star forming regions does have a direct relevance to the 

^ Small jumps may be seen at O.SMq and IMq as the function steps 
suddenly between different modes rather than being smooth, however 
a smooth function merely evens-out these jumps but does not change 
the overall appearance of the IMF. 



origin and form of the IMF. Given the apparent universality of 
the IMF across a wide range of star forming environments (e.g. 
Kroupa 2002) we are presented with two options. Firstly, that 
the mechanism(s) that produce the IMF are fundamentally dif- 
ferent in different environments, but they always produce the 
same outcome. Or, secondly, that there is a single, underlying, 
mechanism that produces the IMF in all environments. The lat- 
ter possibility appeals due to its simplicity, and would suggest 
that the form of the CMF is the driving factor in establishing 
the form of the IMF, and that the form of the CMF is roughly 
the same in diffuse and clustered regions (even if the cores 
themselves are different in spacial size). Indeed, simulations of 
turbulence always seem to produce roughly log-normal CMFs 
whatever the environment. 



6. Conclusions 

We have examined the relationship between the core mass 
function (CMF) and the stellar initial mass function (IMF). We 
use the Orion CMF from Nutter & Ward-Thompson (2007) as 
a 'standard' which we fit using a log-normal distribution. We 
note that this CMF is not dissimilar to the stellar (Kroupa 2002) 
IMF shifted upwards in mass by a factor of 8 (see also Alves 
2007). We randomly sample cores from the CMF and assumed 
that each core produces a certain number of stars with a random 
distribution of masses between the components. 

The canonical IMF is reproduced very well by a scenario in 
which every low-mass cores fragment into binaries, and high- 
mass cores fragment into a multiple system with a ratio of 
binaries-to-triples of 3 : 1 (see e.g. Goodwin & Kroupa 2005) 
and a star formation efficiency (SFE) of ~ 30%. Dynamical 
disruption (Kroupa 1995a,b; Goodwin & Whitworth 2007) of 
systems then evolves the initial binary fraction of unity into the 
field population. 

We find that a scenario in which low-mass stars preferen- 
tially form single systems (e.g. Lada 2006) cannot reproduce 
the observed IMF from a log-normal CMF. Firstly, the slope 
of the high-mass IMF is too steep. Secondly, and most se- 
riously, this model cannot reproduce the correct numbers of 
brown dwarfs to high-mass stars. The best-fit to the canonical 
IMF is found when the SFE is only ~ 15%. Such a low SFE 
is required, as the only way in which brown dwarfs may be 
produced in significant numbers is through the formation of a 
single brown dwarf from a core. Higher SFEs are required to 
produce sufficient high-mass stars, however such SFEs signifi- 
cantly under-produce brown dwarfs and low-mass stars. 

A lingering question is the value of the star formation ef- 
ficiency that must be applied to fit the IMF. The best-fit value 
of e in the fully multiple model suggest that only ~ 30% of the 
mass in a core ends-up in the stars which that core forms (a 
similar value for the SFE is found by Alves et al. 2007). This 
seems a very low value and may suggest that the determina- 
tions of the absolute core masses are wrong. Another possibil- 
ity is that feedback from jets is far more efficient than previ- 
ously thought and manages to disperse most of the gas initially 
in the core. A final possibility is that we are not observing 'typ- 
ical' cores which produce the IMF and that the observed CMFs 
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Fig. 3. Low-mass single star model with A^* -2:1 single:binary for eMc < 0.5Mq, 1 : 1 single:binary for 0.5 < cMc/Mq < 1, 
and 3 : 1 binary:triple for eMc > IMq, with e - 0.15. The initial MSMF (dashed line), and the SSMF (red dotted line) compared 
to the canonical IMF (thin solid lines). 



will produce somewhat top-heavy IMFs (cf. Taurus, Goodwin 
et al. 2004c). 

We conclude that a model in which all stars and brown 
dwarfs form in multiple systems from a log-normal core mass 
distribution provides a very good fit the observed IMF. 
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